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ABSTRACT 

Locating planets in Habitable Zones (HZs) around other stars is a growing field in contem¬ 
porary astronomy. Since a large percentage of all G-M stars in the solar neighborhood are 
expected to be part of binary or multiple stellar systems, investigations of whether habitable 
planets are likely to be discovered in such environments are of prime interest to the scientific 
community. As current exoplanet statistics predicts that the chances are higher to find new 
worlds in systems that are already known to have planets, we examine four known extrasolar 
planetary systems in tight binaries in order to determine their capacity to host additional hab¬ 
itable terrestrial planets. Those systems are Gliese 86, y Cephei, HD 41004 and HD 196885. 
In the case of y Cephei, our results suggest that only the M dwarf companion could host ad¬ 
ditional potentially habitable worlds. Neither could we identify stable, potentially habitable 
regions around HD 196885 A. HD 196885 B can be considered a slightly more promising tar¬ 
get in the search for Earth-twins. Gliese 86 A turned out to be a very good candidate, assuming 
that the system’s history has not been excessively violent. For HD 41004 we have identified 
admissible stable orbits for habitable planets, but those strongly depend on the parameters of 
the system. A more detailed investigation shows that for some initial conditions stable plane¬ 
tary motion is possible in the HZ of HD 41004 A. In spite of the massive companion HD41004 
Bb we found that HD41004 B, too, could host additional habitable worlds. 

Key words: (stars:) binaries: general, planets and satellites: dynamical evolution and stability, 
astrobiology, planets and satellites: individual: HD 41004, Gliese 86, HD 196885 


1 INTRODUCTION 


The past decade has seen a staggering growth in number and 
diversity of known extrasolar planets. The discovery of plane¬ 
tary companions in and aroun d binary stars systems is certainly 
a hig h light in this respec t (e.g. Lathamet ail l 19891: iThorsett et alj 
1993: lDovle etalJ |201ll : IWelsh et alj |2012|: iKostov et all I20l4 


Howard et all 20l4) . We now know of more than 60 multiple star 


systems that harbo r one or more extrasolar planets jSchneider et alj 
I 2 OIII: iReinl l20l4n Some of those systems, e.g. KIC 9632 895 
and Gliese 667 dWelsh et alj2oT^ : lAnglada-Escude et alj2012l) are 
even hosting planets in their respective Habitable Zones (HZ), i.e. 
the regions whe re liquid water can exist on the surface of an Earth¬ 
like planet (e.g. Kasting et al.ll993l : lKaltenegger & Sasselovl201ll : 


iKopparapu et al. 201 4 . There are also works, which showed that 
planet formation is possible in binar y systems, even if the two 
stars h a ve a quite small separation te.g. lHaghighipour & Raymond 
(l2007h . iRafikovI ( I20l4 l. It is interesting to note that nearly one 
third of the currently known exoplanets in S-typ41 configurations 


are part of multi-star multi-planet systems, i.e. of sy stems where 
more than one planet is orbiting the same star (see e.g. lRei nl( l20l4 
or the Binary Catalogue of ExoplanetB)- Current statistics of the 
Kepler candidates exoplanet population furthermore suggest that 
about 4 6% of all planets discovered so far reside in multiplanet 
systems dBurke et alj20l4 . Despite the fact that the planetary mul- 
tiplicity fraction in multiple-star systems is not well constra ined 
dThorp et alJl20l4lGinski et'^l20l4lArmstrong et alj|20l4 . the 
high ratio of multiplanet systems in Kepler data may foster hopes to 
discover new planet candidates in multiple star associations that are 
already known to host a giant planet. Such additional planets would 
have to be small or far out from their host star on long periodic or¬ 
bits, o therwise they would have been detected by now dEggl et alj 
l2013h . In this study we investigate whether it is possible for binary 
star systems with a known gas giant (GG) to host additional terres¬ 
trial planets in their respective HZs, since the latter are attractive 
observational targets. In particular we selected the following four 
binary systems with separation of approximately 20 au, since the 
strongest gravitational perturbations can be expected in such sys¬ 
tems. 


* E-mail: funk@astro.univie.ac.at; 

' http://www.openexoplanetcatalogue.com/, http://www.exoplanet.eu/ 

^ In this configuration the planet moves in an orbit around one of the stars. 


http://www.univie.ac.at/adg/schwarz/multiple.html 
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• 7 Cephei: The 7 Cephei binary system consi sts of a K1 III-IV 
star an d a dwarf companion of spectral type M4. iNeuhauser et alj 
ll2Q07h determined the orbital parameters of the binary system 
(as shown in Table [p. Collecting RV measurements spanning 


many years 


J Col l 
2011 

2011 


EndletalJ i 2011 ) updated the orbital parameters and 


iReffert & Ouirrenbachl ( 2011 ) used Hipparcos data to search for 
the astrometric signature of planets and brown dwarfs. For the sys¬ 
tem 7 Cephei the signature of the substellar companion was strong 
enough to give reasonable constraints on inclination, ascending 
node and mass. They claim that 7 Cephei b is most likely a brown 

dwarf. _ 

• Gliese 86: lOueloz et al.l ( I 2 OOOI) discovered a 4 Mi„„ plane t 
around the K1 V star Gliese 86. One year later IeIs et alJ l l200lh 
announced a substellar companion with a mass of approximately 
50 Mj ,i„ in a distance of 18.75 au. Later iMugrauer & Neuhau^ 
( I 2 OO 5 ) confirmed that the secondary is a White Dwarf (WD) and 


not a brown dwarf. Wit h the help of new observations using the 
Hubble Space Telescope lFarihi et al.l ( 1201 3h could proof that Gliese 
86 B is a WD and could give modified data for the system. 

• HD 196885: The system HD 196885 consists of an F8 V 
star with 1.3 a Ml V companion and a planet orbitin g 

the primary (see e.g . lOiauyin eLaj 12006). Ichauvin et al. |2007h, 


tne primary (see e.g . |Unauvm_et_aU (), |Uiauvm_et_alJ^lzyt^, 
ICorreia et ^20081) or Fischer et alJ ' 2001^ )). Chauvin et al.l ^201111 


considered astrometric data as well as all data from past RV sur¬ 
veys (CORAVEL, ELODIE, CORALIE and Lick, over more than 
26 years) to constrain the physical and orbital properties of the HD 
196885 system and could find a unique best solution, w hich is sum¬ 
marized in Table[T] Furthermore lChauvin et al.l (1201 ll) investigated 
the stability of the system and showed that the system is more stable 
in a high mutual inclination configuration. 

• HD 41004: The system HD 41004 comprises a K1 V primary 


of 20 au. ISantos et al. 

d2002r) discovered a brown dwarf orbiting 

HD 41004 B and later 

Zucker et alJ (|2004|) found a GG around HD 


41004 A. 

In the following we investigate, whether dynamical stability allows 
for an additional terrestrial planet in the HZs of the four considered 
systems. 


2 THE HABITABLE ZONES 

The definition of the HZ is a highly interdisciplinary topic. Im¬ 
portant parameters are not only the mass and atm osphere of the 
potentially habitable planet iKopparapu et alj2014h , but also prop¬ 
erties of the host star like activity and radiation (e.g. dLammer et al.l 
l2009t) ) and last but not least the dynamical stability of all involved 
bodies play a significant role. The evolution of orbital parameters 
can also be a determining factor. Large variations in a planet’s or¬ 
bital eccentrici ty, for instance, might chan ge the conditions for hab¬ 
itability (e.g. dWilliams & Pollaral2002h ). Considerations of dy¬ 
namical stability are of special interest for multiple star planet 
systems due t o potentially larg e gravitational per turbations of the 
planet’s orbit dlaime et alj20 A . lEggl et alJd2012h investigated the 
influence of stellar companions with different spectral types on the 
insolation a terrestrial planet receives orbiting a Sun-like primary. 
They used a time-independent analytical estimate to calculate the 
borders of the HZ and found out that there is a strong dependence of 
permanent habitability on the binary’s eccentricity. With four bod¬ 
ies in a system, special attention has to be given to the assessment 
of habitability. The following section contains a brief discussion on 
the related issues. 


2.1 Calculating the borders of the HZs 


In this section we will shortly summarize how to determine the HZ 
for the four investigated binary star systems. Several approaches 
can be found in literature how to calculate the borders of th e 
circumstellar HZ in binary star sy stems (see e.g. ICuntzl d2014h . 
iKaltenegger & HaghiehipouJ d2013h ). We choose the method pre¬ 
sented in Eggl et al] ( ITOl^. because it com bines the radiative 
and the dynamical effects. lEggl et alj (1201 2h investigated differ¬ 
ent binary-planet configurations and examined the implications of 
stellar companions with different spectral types on the insolation 
a terrestrial planet receives orbiting a Sun-like primary. The time- 
indepe ndent analytical estimates that were presented in lEggl et al.l 
( 1201 2[) are based on three types of HZs: 


• Permanently Habitable Zone (PHZ): The PHZ is the region 
where a planet always stays within the insolation limits of the cor¬ 
responding HZ. 

• Extended Habitable Zone (EHZ): In contrast to the PHZ 
parts of the planetary orbit lie outside the HZ. Yet, the binary planet 
configuration is still considered to be habitable when most of its or¬ 
bit remains inside the HZ boundaries. 

• Averaged Habitable Zo ne (AHZ): Following the argument 
bv lWilliams & PollardI ( |2002|) . this category encompasses all con¬ 
figurations which allow for the planet’s time-averaged effective in¬ 
solation to be within the limits of the HZ. 


lEggl et al] ( l2012h found that the AHZ is almost independent of the 
binary’s ecce n tricity and coincides well with the HZ defined by 
iKasting et alj ( Il993ri . while the PHZ and EHZ shrink drastically 
with higher binary eccentricities. The gravitational influence of the 
second star injects eccentricity into the planet’s orbit, which in turn 
leads to considerable changes in planetary insolation due to closer 
encounters with the host-star. 

The HZs of 7 Cephei, HD 41004, HD 196885 and Glie se 86 were 
determ ined using the analytical method present ed in lEggl et al 


2OI2I) with updated effective insolation values dKopparapu et al 


120146 . The results are shown in Figures[T] 


Blue regions denote the extent of the PHZ, green regions cor¬ 
respond to EHZs and yellow represents AHZs. Red regions are not 
habitable, but dynamically stable, while the shaded regions are dy¬ 
namically unstable. The black lines corresponds to the classical HZ 
where the second star is completely neglected. The initial condi¬ 
tions for the binaries are summarized in Table [T] Due to the large 
orbital periods of S-type binary stars, the eccentricity of the binary 
systems is often poorly constrained. Hence, results for three differ¬ 
ent eccentricity values of the binary’s orbit (eB,„ = 0.2, 0.4 and 0.6) 
are provided. 

The borders of the PHZ, EHZ and AHZ are summarized in Ta¬ 
ble |2] For all four systems one can see that the PHZ shrinks 
with h igher eccentricities. A comparison of our work with 
that of IKaltenegger & Haghighipouj ( 1201 3h can be done for the 
system HD 196885, which was investigated in both studies. 
IKaltenegger & Haghighipo^ (1201 3l) determined the following bor¬ 
ders for the HZs of components A and E0: HZa '. 1.454 au - 2.477 


Considering the maximum flux of the companion, which corresponds 
to the closest approach of the stars and therefore, to the strongest in- 
fluence on the HZ. Since no dyna mical perturbations are considered in 
IKaltenegger & Haghighino^ j2013h the HZ borders for the minimum and 
maximum flux values are almost the same (deviations between 0.002 au and 
0.008 au). 
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Table 1. Keplerian orbital parameters and errors (if given) of the four known tight binary systems, hosting an extrasolar planet. The orbital elements of 
component B are given with respect to component A and the elements of the planets are given with respect to the respective central star. Given are the Name, 
semi-major axis (a), eccentricity (e), argument of perihelion (tu). Longitude of the ascending node (SI), mass, spectral type and age. 


Name 

a [au] 

e 

0 ) [deg] 

SI [deg] 

mass 

Sp.Type 

Age [Gyr] 

Ref. 

y Cephei A 

- 

- 

- 

- 

1.4M5„„ 

±0.12 

K1III-IV 

6.6 

Neuhauser et al. 120071 

y Cephei B 

20.18 

0.4112 

162.0 

18.04 

0.409 Ms„„ 

M4 


Neuhauser et al. (20071 


± 0.66 

± 0.0063 

±0.40 

±0.98 

±0.018 




y Cephei Ab 

2.05 

0.049 

94.6 

37.5 

1.85 Mjup 

- 


Endletal. 120111. 

Reffert & Ouirrenbach (201 1) 


±0.06 

± 0.034 

±34.6 


±0.06 




Gliese 86 A 

- 

- 

- 

- 

0.8 Msun 

KOV 

2.03 

Farihi 6^^200) 

Gliese 86 B 

27.8 - 69.8 

0.0-0.61 

- 

- 

0.59 Msu„ 
±0.01 

DQ6 


Farihi et al. 12013) 

Gliese 86 Ab 

0.11 

0.046 

270.0 

- 

4.0 Mjup 

- 


Oueloz et al. (20001 



± 0.004 

±4 






HD 196885 A 

- 

- 

- 

- 

1.3 Msun 

F8 V 

2.0 

Chauvh^^^20^) 








±0.5 


HD 196885 B 

21 

0.42 

-118.1 

79.8 

0.45 Msun 

Ml V 


Chauvin et al. (20^) 


±0.86 

±0.03 

±3.1 

±0.1 

±0.1 




HD 196885 Ab 

2.6 

0.48 

93.2 

- 

2.98 Mjup 

- 


Chauvin et al. (20^) 


±0.1 

±0.02 

±3.0 


±0.05 




HD 41004 A 

- 

- 

- 

- 

0.7 Msun 

K1 V 

1.64 

Santos et al. (2002) 

HD 41004 B 

20 

0.4 

- 

- 

0.42 Msun 

M2 V 

1.56 

Roell et al. 120121. Zucker et al. 120041 








±0.8 


HD 41004 Ab 

1.6 

0.39 

97 

- 

2.54 Mjup 

- 


Roell et al. 120121. Zucker et al. 120041 



±0.17 

±31 


±0.74 




HD 41004 Bb 

0.0177 

0.081 

178.5 

- 

18.37 My„p 

- 


Zucker et al. 12004), Santos et al. 12002) 



±0.012 

±7.8 


±0.22 





GammaCeph S-Type A 



0 1 2 3 4 5 6 7 

planet semimajor axis (ap) 


HD41004 S-Type A 



0 0.6 1 1.5 2 2.5 3 3.5 4 

planet semimajor axis (ap) 


HD196885 S-Type A 



0 0.5 1 1.5 2 2.5 3 3.5 4 

planet semimajor axis (a ) 


Gliese86 S-Type A 



0 0.6 1 1.5 2 2.5 3 3.5 4 

planet semimajor axis (ap) 


Figure 1. Circumprimary HZs for the systems (from left to right) y Ceph, HD 41004 , HD 196885 and Gliese 86, where blue gives the PHZ, green the EHZ 
and yellow the AHZ. Red regions are not habitable, but dynamically stable. Shaded regions are dynamically unstable. The black lines gives the extent of the 
HZ, when ignoring the secondary star. See online version for color figures. 


au; HZb'. 0.260 au - 0.491 au. Comp aring our results, we find that 
the PH Z lies within the borders of iKaltenegger & Haghighipouj 
( l2013h . which is due to the consideration of the perturbing influ¬ 
ence of the secondary. The EHZ allows eccentric planetary orbits, 
that may le ave the HZ for a short time, which lea ds to quite similar 
results as in iKaltenegger & Haghighipouil j2013h . 

Please note that the analytical method used to determine the pre¬ 
sented HZ borders assumes that the luminosity of the stars does not 
evolve significantly. For stars on the main sequence this is true for 


billions of years. For post main sequence stars, the timescale on 
which those estimates remain valid depends on the speed at which 
the stellar luminosity evolves. Furthermore, the analytic solutions 
are based on eccentricity estimates that do not consider the known 
GG. They merely serve as initial guidelines to determine search 
spaces for additional, potentially habitable terrestrial planets. If the 
GG influences the eccentricity evolution of additional planets sig¬ 
nificantly, a separate assessment of their habitability is required, see 
section[6j2l 
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GammaCeph S-Type B 



0 0.2 0.4 0.6 0.8 1 

planet semimajor axis (ap) 


HD41004 S-Type B 



0 0.2 0.4 0.6 0.8 1 

planet semimajor axis (ap) 


HD196885 S-Type B 



0 0.2 0.4 0.6 0.8 1 

planet semimajor axis (ap) 


Gliese86 S-Type B 



0 0.02 0.04 0.06 0.08 0.1 

planet semimajor axis (ap) 


Figure 2. Same as Figrue[T] only for HZs around the less massive component (B) of the binary star. From left to right: y Ceph, HD 41004 , HD 196885 and 
Gliese 86. Note the difference in the scale of the planetary semimajor axis for Gliese 86. See online version for color figures. 


Table 2. Values of the AHZ, EHZ and PHZ (in au) for the systems y Cephei, HD 41004, HD 196885 and Gliese 86. The HZ borders are provided for the 
primary (A) and the secondary star (B). Since the eccentricity of some of the binaries is not well constrained we give three values for each system corresponding 
to different binary star eccentricities: esm = 0.2, egm = 0.4, and esm = 0.6. 


Name 

inner AHZ 

inner EHZ 

inner PHZ 

outer PHZ 

outer EHZ 

outer AHZ 

y Cephei A {egm = 0-2) 

4.17 

4.36 

4.79 

5.70 

5.70 

5.70 

y Cephei A (esm = 0.4) 

- 

- 

- 

- 

- 

- 

y Cephei A {esm = 0.6) 

- 

- 

- 

- 

- 

- 

y Cephei B {esm = 0.2) 

0.183 

0.183 

0.187 

0.363 

0.377 

0.379 

y Cephei B (esin = 0.4) 

0.185 

0.185 

0.193 

0.355 

0.385 

0.389 

y Cephei B {esi,, = 0.6) 

0.189 

0.191 

0.215 

0.345 

0.421 

0.431 

HD 41004 A {esin = 0.2) 

0.77 

0.77 

0.78 

1.34 

1.36 

1.38 

HD 41004 A {eBm = 0.4) 

0.77 

0.78 

0.80 

1.29 

1.34 

1.39 

HD 41004 A (eg,,. = 0.6) 

0.77 

0.80 

0.85 

1.22 

1.30 

1.40 

HD 41004 B (eg,,, = 0.2) 

0.374 

0.376 

0.378 

0.712 

0.720 

0.726 

HD 41004 B (eg,„ = 0.4) 

0.374 

0.378 

0.382 

0.696 

0.712 

0.726 

HD 41004 B (eg,„ = 0.6) 

0.374 

0.382 

0.392 

0.674 

0.704 

0.730 

HD 196885 A (eB,„ = 0.2) 

1.43 

1.45 

1.48 

2.35 

2.42 

2.48 

HD 196885 A (eein = 0.4) 

1.43 

1.48 

1.56 

2.22 

2.35 

2.50 

HD 196885 A (esin = 0.6) 

1.44 

1.55 

1.74 

2.04 

2.27 

2.39 

HD 196885 B (eem = 0.2) 

0.255 

0.257 

0.257 

0.489 

0.493 

0.497 

HD 196885 B (f>B,„ = 0.4) 

0.255 

0.257 

0.261 

0.481 

0.491 

0.497 

HD 196885 B {egin = 0.6) 

0.257 

0.259 

0.267 

0.471 

0.489 

0.501 

Gliese 86 A (eg,„ = 0.2) 

0.68 

0.69 

0.69 

1.20 

1.21 

1.22 

Gliese 86 A (eg,,, = 0.4) 

0.68 

0.69 

0.70 

1.17 

1.19 

1.22 

Gliese 86 A (eg,„ = 0.6) 

0.68 

0.70 

0.72 

1.12 

1.17 

1.23 

Gliese 86 B {esm = 0.2) 

0.023 

0.023 

0.023 

0.037 

0.037 

0.037 

Gliese 86 B (eg,„ = 0.4) 

0.023 

0.023 

0.023 

0.037 

0.037 

0.037 

Gliese 86 B (eg,„ = 0.6) 

0.023 

0.023 

0.023 

0.037 

0.037 

0.037 


2.2 Stability of the HZs 

In order to ensure the dynamical stability of an additional terrestrial 
planet in the HZ of the four investigated systems we studied the 
dynamical evolution of test-planets within the borders of the 
previous calculated region of the HZ (see Table|3 for each system. 
As a model we used the elliptic restricted four-body problem, 
consisting of the two stars, the known GG and mass-less TP’s, 
which is a good approximation for terrestrial planets. In the case 


of HD 41004 the masses of component B and of the nearby brown 
dwarf were added and put to there common barycenter, since 
test computations showed that the brown dwarf has no perturbing 
influence on the region of the HZ (see section 6.3). In all cases, 
the motion of the planets and the binary star were treated self 
consistently in the framework of Newtonian mechanics. All the 
celestial bodies involved were regarded as point masses. Test 
computations have shown that a time of 10® years is sufficient to 
capture the most relevant dynamical effects. For all integrations 
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Figure 3. Maximum eccentricities for three different eccentricities of Gliese 
86 B (green: eg,,, = 0.0, blue: eg,,, = 0.7, red: egm = 0.8). On the jr-axis we 
show the semi-major axis, covering the HZ. The v-axis gives the maximum 
eccentricity values. See online version for color figures. 

we used the Lie-Series Integration Method (H anslrne ier & Dvorak! 
( ll984h . [Lichteneg^ jl984f) . [Eggl & Dvorak j2010h ) which is also 
capable of dealing with large eccentricities and close encounters 
between bodies. To analyze the data we used the maximum 
eccentricity (ME) that the test planets (TP) orbit reached during the 
integration time. If the TP’s orbit becomes parabolic or hyperbolic 
(ME >1) with respect to it’s central star the orbit is considered to 
be unstable. A further advantage of using ME for the investigations 
of TP’s in the HZ is that information on the extent of the PHZ 
becomes directly available. 

The orbital parameters of the investigated four systems are 
summarized in Table [T] The mass-less TP’s were started equally 
distributed (Art = 0.01 au) within the borders of the HZ. All other 
orbital elements of the test planets (e, i, ai, fl and M) were initially 
set to zero. 

3 7 CEPHEI 

The 7 Cephei system has already been investigated in detail con- 
cerning the dynami cal stability in its c ircum primary HZ (see e.g. 
[Dvorak et al.l ( l2003l) or lHaghighinoud ( l2006h ). However, the con¬ 
tribution of the secondary to the extent of the HZ w as completely 
neglec ted in those studies. The numerical study by [Dvorak et al.l 
( [2OO3I) examined the area between 0.5 and 1.85 au concerning sta¬ 
bility. The stability maps showed a stable region between 0.5 and 
1.2 au and a chaotic area for lager semi-major axes of the test- 
planet. Within the stable zone an arched chaotic band was found 
which results from the combi ned perturbation of th e giant planet 
and the secondary star (see e.g. lPilat-LohingeJ 1 I2OO5I) ). The authors 
studied also the area around 1 au in detail using several masses (up 
to 90 Earth masses) for the small planet moving in this area. Even 
if they talked about habitability when discussing the results it is 
well known that the studied region does not corresp ond to the HZ 
of 7 C ephei. This was also criticized in the study bv [HaghighiDouil 
( I2OO6I) where the area between 0.3 and 4 au was studied concerning 
long-term stability and habitability. Taking into account the stellar 
type Haghighipour defined the HZ as the area between 3.05 and 3.6 
au for which he could not get long-term stable planetary orbits. 
Updates for the HZ around 7 Cephei A show that the habitable 
limits are shifted beyond the orbital stability limits of the current 
system (cb,,, = 0.4112). Only stable orbits in the HZ around the sec¬ 


ondary are possible. As no planet is known to orbit 7 Cephei B the 
HZ borders given in Table 1 can serve as an observational guideline 
provided one assumes that stable circulation pattern s have formed 
on the additional terrestrial planet jWang et al.l2014[) and tidal heat¬ 
ing remains negligible. 


4 HD 196885 

For the system HD 196885 our results provide constraints on the 
possible eccentricity of the binary. Since we know a GG at 2.6 au 
the given eccentricity of the binary (e = 0.42) is quite close to the 
system’s stability limit. Comparing Tables[T]and[2]one can see that 
the known planet in the systems HD 196885 and HD 41004 orbits 
close to or even inside the HZ. For HD 196885 the GG is located in 
the middle of the HZ. Given its high eccentricity of e = 0.48 it will 
cause dynamical instability for any planet in this region. Test com¬ 
putations showed that even if the known planet had no eccentricity 
(e = 0) no additional stable orbits would be possible. Similar to 
7 Cephei, the only possibility for HD 196885 to host additional po¬ 
tentially habitable planets is to have them orbit the M-dwarf com¬ 
panion. 


5 GLIESE 86 

A detailed dynamical study on Gliese 86 has already been per¬ 
formed bv IPilat-Lohinger & Funid ( |2006|) . They found that the GG 
does not dynamically influence the region of the HZ. Since the or¬ 
bital parameters of Gliese 86 as well as the method of determining 
the HZ have changed we shall check their results for robustness. To 
this end, we compared the analytically obtained HZ (without GG, 
Figure [TJ with numerical integrations of the whole system (includ¬ 
ing the GG). Since the binary star’s orbit is not well constrained, we 
integrated this system for different eccentricities. Figure [3 shows 
the results for e = 0.0, 0.7 and 0.8. The jc-axis the semi-major axis, 
covering the best-case HZ (eB,„ = 0, see Table|2]l. The y-axis gives 
the maximum eccentricity value. It turns out that for binary eccen¬ 
tricities up to e = 0.7 the TP’s maximum eccentricity reaches values 
up to 0.2 only. Thus we can conclude that stable habitable planets 
are possible in the Gliese 86 system even if the binary star’s or¬ 
bit is highly elliptic. Only for binary eccentricities higher than e = 
0.7 the secondary’s influence on the test planets in the region of 
the HZ becomes critical. Comparing these results with the analyt¬ 
ically obtained HZ borders (see Table [2] and Figure [3 one can see 
that similar results can be obtained for the test planets eccentric¬ 
ity, which means that the known GG has no significant influence 
on additional planets in the HZ. Therefore, the analytically derived 
borders of the HZ remain valid even for the whole system of Gliese 
86. 

Nevertheless, we are faced with two problems when talking 
about habitable planets in Gliese 86. The first issue is related to the 
fact that Gliese 86 B is a WD. This suggests a cataclysmic phase 
in the evolution of the Gliese 86 system with strong variations in 
stellar insolation and mass loss. Whether an Earth-like planet could 
retain any sort of habitability during such events or can be capture d 
afterward is discussed in e.g. [Agoll ( l201lh . [Bames & Hell^ ( [2013h . 
[Veras & Gansickel l l2014l) . The second issue is the small distance 
of the known GG to its host star. This could imply that migration 
processes occurred either during the system’s formation, or during 
the cataclysmic phase. If the GG formed beyond the Gliese 86 A’s 
snowline and migrated inward, it had to cross the HZ. This would 
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have destabilized any possible terrestrial planet in the circumpri- 
mary HZ. 

However, two scenarios are thinkable that would still allow 
for the existence of potentially habitable worlds around Gliese 
86 A: a) the planet w as built after the migration of the GG (e.g. 
[Raymond et al.|[200^ . or b) the planet formed outside of the GG, 
was trapped in a mean motion resonance and migrated inward to¬ 
gether with the GG. Such a scenario would allow for a large pri¬ 
mordial water content of the rocky planet. 

While it is clear that Gliese86 B has gone through a cata¬ 
clysmic phase, it is intriguing to note that a circumstellar habitable 
planet around the B component seems possible. Although a WD, 
Gliese 86 B does have a circumstellar region that would allow for 
liquid water on an Earth-like planet’s surface. This region is very 
small and very close to the star, as can be seen from Figure[2 How¬ 
ever, tidal forces do not yet outweigh the self gravity of Earth-like 
planets at this distance, since the fraction of accelerations due to 
the planet’s self gravity asa and tides aj acting on a particle at the 
top of the atmosphere 100km above the Earth’s surface (Karman 
line) is still much larger than unity: 


— ^ ~ 40 

flr 


( 1 ) 


for an Earthlike planet orbiting Gliese86 (B) in the HZ. Here, nip 
and Rp denote the planet’s mass and radius to the top of the atmo¬ 
sphere, respectively. Furthermore, m* is the stellar mass and r the 
planet’s orbital distance from the host star. Hence, if an Earth-like 
planet had survived the cataclysmic event, or if it had formed later 
on, it could retain its atmosphere. If it orbits in the HZ it could 
be potentially habitable - from the point of view of insolation. A 
tidally locked rotation state and non negligible tidal heating would 
be likely in such environments. Such factors require a more detailed 
analysis that lies beyond the scope of this article. 


6 HD 41004 

6.1 Dynamical Study 


HD41004 is a very interesting system, because the known giant 
planet’s orbit is close enough to the HZ to significantly influence 
additional potentially habitable planets. First investigations con¬ 
cerning the stability in t he reg ion of the HZ of HD41004A were 
done by IPilat-Lohingej ( 12005 ) showing mean motion and secular 
perturbations due to the interactions of the discovered GG and the 
secondary star. The study shows different dynamical maps for var¬ 
ious semi-major axes of the gas planet where the appearance of 
mean motion and secular perturbations strongly depend on the po¬ 
sition of the GG. In some cases these perturbations interfere with 
the HZ. Due to uncertainties and changes in the orbital parameters 
new computations of this pla netary system were performed (see 
IPilat-Lohinger & FunkI ( 120 lOh ) for the different planetary configu¬ 
rations. This study revealed the possible maximum eccentricity of 
the binary stars for the different orbital parameters given by the ob¬ 
servations. 

Even if preliminary calculations of additional hypothetical terres¬ 
trial planets in this binary system suggested a zero probability for 
other planets in the HZ due to stability issues slight variations in 
the nominal ec centricity and semi-major a xis change this picture 
completely (see lKlat-Lohinger et alj ( l2012li ). 

The nominal eccentricity of the (GG) is e = 0.39 ± 0.17 the planet’s 
orbit is rather badly constrained. The same holds for the binary 


star’s orbit. Therefore, we investigated whether any combination of 
permissible eccentricities of the system could yield stable zones in 
the HZ. Fixing the remaining orbital parameters of HD 41004 as 
given in Table[T] we studied this system with decreased eccentrici¬ 
ties of both the secondary (eg,,, = 0.2) and the GG (eoc = 0.2). 

For all integrations we started mass-less test planets in a region be¬ 
tween 0.2 and 1.3 au, which covers most of the circumprimary HZ 
of HD 41004, see Table Since the region of the HZ is unsta¬ 
ble when using the published orbital parameters of HD41004 we 
continued the computations for smaller semi-major axes in order to 
find the border of stability. We stopped our integrations at 1.3 au 
due to the presence of a GG at 1.64 au. The results are summarized 
in Figure [4] where on the x-axis the semi-major axis of the test 
planets is shown and on the y-axis their inclination with respect to 
the binaries plane, which could also be an important factor. As one 
can see in the top panel of Figure [4] upper graph no stable motion 
is possible in the region of the HZ for the given orbital parameters. 
When the secondary’s eccentricity is decreased (e^,,, = 0.2, Figure 
[4] middle graph) the border of stable motion shifts just slightly out¬ 
wards, but here one can also see the influence of the secondary, 
caused by a secular resonance between a x 0.25 and 0.5 au. This 
secular resonance becomes much bet ter visible in Figure [4 lowe r 
graph and is investigated in detail by IPilat-Lohinger et ( 2014h . 
In Figure [4] lower graph one can see clearly that the dominant in¬ 
fluence on the HZ comes from the GG, since the stable region is 
shifted significantly outwards for a lower eccentricity of the GG 
(fiGG = 0.2). Nevertheless the region of the HZ is still - apart from 
a few orbits just on the inner edge - mainly unstable. 

Thus, we investigated the influence of changes in the semi-major 
axis of the known GG. Since the region of the HZ turned out to 
be most stable at the inner edge of the HZ we fixed in a first step 
the semi-major axis of the test planet at 0.8 au and varied the semi¬ 
major axis of the GG between 1.5 and 5 au. The integrations were 
done for 5 x 10® years for eg,,, = 0.2 and 0.4 and for ecc = 0.2 
and 0.4. This test integrations showed, that long term stable motion 
in the HZ is just possible if eB,„ = Ccq = 0.2. For higher eccen¬ 
tricities no stable motion of the test planets could be found - if the 
GG’s semi-major axis is low it perturbs the test planet and when 
the semi-major axis of the GG becomes higher it is perturbed it¬ 
self by the secondary, which again leads to perturbations on the test 
planet. For the best case scenario (eB,„ = Sgc = 0.2) we constructed 
a grid of initial conditions, where on the x-axis the semi-major axis 
of the test planets and on the y-axis the semi-major axis of the GG 
are shown. The results are displayed in Figure where the white 
line indicates the current position of HD 41004 b. Here some sta¬ 
ble regions can be seen also for the given semi-major axis of HD 
41004 A b. A detailed investigation of the dynamical stability i n 
the system HD 41004 can be found in IPilat-Lohinger etd] ( l2014h . 
but here we will concentrate on habitability considerations. 


6.2 The HZ of HD 41004 A 

Following the results of the previous section, we aim to ascertain 
that the stable regions in Figure |5] also fulfill the criteria for habit¬ 
ability. In section im we defined the three body HZ including the 
influence of the secondary but neglecting the G G. Hence, we have 
to take the influence of the GG into account now. lEggl et alj ( l2012h 
already showed that the PHZ of a specific planet in a binary sys¬ 
tem is a function of the pericenter and apocenter distances from the 
primary star. In hierarchical systems, those are mostly determined 
by the planet’s eccentricity evolution. In our case the eccentricity 
of the test planet does not only change because of the influence of 
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Figure 4. Semi-major axis of the test planets versus their inclination. Upper 
graph: given orbital parameters for the system HD 41004; middle graph: 
like upper graph, but for esm = 0.2; lower graph; like middle graph, but 
for eQQ = 0.2. The color code coiTesponds to the maximum eccentricity of 
the test planets, where violet shows unstable motion. The white vertical line 
shows the inner border of the PHZ. See online version for color figures. 




the secondary, but also due to the perturbations of the GG. Using 
the output of the dynamical investigations we define the region of 
PHZ via the maximum eccentricity values. Results are presented in 
Figure [h] which shows graphs similar to Figure From Figure 
it becomes clear that stable configurations for a GG’s semi-major 
axis of 1.64 au are not in the PHZ. Nevertheless, we could find 
some stable, habitable configurations, e.g. when the GG is at 2 au 
with an eccentricity of 0.2. 


Figure 5. Semi-major axis of the test planets versus semi-major axis of the 
gas giant for an eccentricity of the planet of bcq = 0.2 and an eccentricity 
of the binary of esin = 0.2. All other orbital elements are as shown in Ta- 
ble[T] the test planets move on a planar orbit. The color code corresponds to 
the maximum eccentricity of the test planets, where violet shows unstable 
motion. The white line indicates the current position of HD 41004 A b. See 
online version for color figures. 
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Figure 6. Same as Figure|^ but the color code shows permanent habitability 
(PHZ, blue), and non-habitable regions (red). The latter are dynamically 
stable but not in the PHZ. Gray areas denote dynamically unstable regions. 
See online version for color figures. 



Finally, we investigated the influence of the binaries mutual dis¬ 
tance. To this end, we fixed the test planets position at 1 au and 
varied the semi-major axis of the GG and of the secondary. The 
test computations were again done for different eccentricities of the 
GG and the secondary (bqq = 0.2 and 0.4; esm = 0.2 and 0.4). All 
other initial conditions are provided in Table [T] The corresponding 
results are shown in Figure |7] where the color code again corre¬ 
sponds to the potential habitability. For a GG’s eccentricity of 0.2 
(left graphs) stable habitable configurations are possible for a GG’s 
distance larger than » 2 au. Interesting is that the GG’s distance 
needs to be larger for larger distances of the secondary to allow 
stable configurations in the PHZ (ace ~ 2 au for agm = 10 au 
up to Occ ~ 3.8 au for agj„ = 40 au). This is caused by the per- 
turbing secular frequency of the GG (g = gcc^ g c™ be found in 
IPilat-Lohinger et ^ ( 1200 Sh ). which is visible as the gray diagonal 
region in Figure |7] (also visible as the arc like structure in Figure |4] 
bottom g raph). A detailed investigati on of this phenomenon can be 
found in IPilat-Lohinger et ^ ( l2014h . However two mean motion 
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resonances between the GG and test planet allow stable habitable 
configurations also within the secular resonance of the secondary, 
namely the 3:1 mean motion resonance at ace = 2.08 au and the 
4:1 mean motion resonance at 2.52 au. 

For a GG’s eccentricity of 0.4 (right graphs) stable, habitable con¬ 
figurations are possible for a GG’s distance of more than » 3.3 au. 
Once more, the perturbing influence of the aforementioned secular 
resonance becomes evident. 


6.3 The HZ of HD 41004 B 

HD 41004 B is itself a binary system with a GG (HD 41004 Bb) 
orbiting its host star at a distance of only roughly 0.02 au, see Ta¬ 
ble [T] Since the three body HZ (star star planet) for the B compo¬ 
nent encompasses the region between 0.347-0.730 au, the question 
arises whether HD 41004 Bb can significantly influence an Earth¬ 
like planet and change the borders of the analytically determined 
HZ. To answer this question, we numerically calculated the varia¬ 
tion in orbital elements of a massive Earth-like planet on an initially 
circular orbit at the inner edge of the HZ of HD 41004 B. Assum¬ 
ing that this system is coplanar, the ME of the Earth-like planet at 
the inner edge of the HZ only grows up to ~ 0.011, and the 
variations in its semimajor axis remain below 10“^ au. Hence, we 
conclude that the GG HD 41004 Bb does not influence planets in 
the HZ significantly and the HZ borders remain accurate. This in 
turn means that it is not impossible to have potentially habitable 
planets in stable orbits around HD 41004 B. 


7 SUMMARY AND DISCUSSION 

The aim of our work was to investigate the possible existence of 
additional terrestrial planets in the Habitable Zones of four known 
close binary systems known to host giant planets. The selected sys¬ 
tems were y Cephei, Gliese 86, HD 41004, and HD 196885. In 
a first step we calculated three body HZs for all systems excluding 
the influence of the known GG. Our investigation of circumprimary 
HZs showed that in two of the systems (HD 41004 and HD 196885) 
the GG are situated within or near the three body HZs, while in the 
third system (Gliese 86) the distance between the circumprimary 
HZ and the GG is large enough that the influence of the latter be¬ 
comes negligible. The current orbital configuration of the y Cephei 
system does not allow for stable planetary motion in the circumpri¬ 
mary HZ. 

Numerical investigations for the system Gliese 86 showed that 
planets in the Habitable Zone are not so much atfected by the 
known GG. The effect of the orbital eccentricity of the secondary 
star is more important, but only if > 0.7. Regarding the sys¬ 
tems HD 41004 and HD 196885 we could show that the system 
HD 196885 does not allow any stable orbit in its HZ, even not if 
the eccentricities of the GG and the binary would be smaller. 

The situation is better for the system HD 41004. Although, first 
integrations showed also that no stable orbits are possible in the re¬ 
gion of the HZ, here small changes (within the given error bars) in 
eccentricity and semi-major axis of the GG and the binary could 
lead to at least some small stable regions in the HZ. Therefore we 
can conclude, that HD 41004 although not very good, is still the 
best candidate for stable motion in its HZ among the four known 
close binary systems hosting a extrasolar planet. In all cases, HZs 
around the less massive component of the double star seem to be 
capable of hosting dynamically stable terrestrial planets. 
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The test-planet was always at 1.0 au; all other initial conditions were set to the values given in Table[T] The color code shows the habitability, where blue 
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figures. 
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